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Q , We compute the effects of the centrifugal pressure supported shock waves on the emitted 

' spectrum from an accretion disk primarily consisting of low angular momentum matter. 

' Electrons are very efficiently accelerated by the accretion shock and acquire power-law 

^ ^ distribution. The accelerated particles in turn emit synchrotron radiation in presence 

of a stochastic magnetic field in equipartition with the gas. Efficient cooling of the 
electrons by these soft photons reduces its temperature in comparison to the protons. We 
explore the nature of the broadband spectra by using Comptonization, bremsstrahlung 
rN( and synchrotron emission. We then show that there could be two crossing points in a 

^ , broadband spectrum, one near ~ lOkeV and the other ~ 300 — 400KeV. 
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1. Introduction 

It is more than fifty years since the study of the accretion fiows on to gravitat- 
ing compact objects began. Bondi^ first showed that spherically symmetric matter 
would pass through a sonic sphere before falling on a sufficiently compact object. 
However, the matter has a large radial motion. Close to the black hole, velocity 
is very high, and thus the density is low. As a result, the flow has a very low ra- 
diation efficiency. Discovery of quasars and active galaxies in this decade required 
that the efficiency be improved. For this, several important steps were taken. First, 
Shvartsman^ introduced dissipation due to entangled magnetic field. Shapiro'^"* 
computed the degree of dissipation and found that the luminosity is increased by a 
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significant amount. However, this was still not sufficient to explain quasar luminos- 
ity. Shakura-Sunyacv^ introduced Keplerian thin disk models of the active galaxies. 
These are efficient and did indeed explain the 'big blue bump' in UV/EUV region 
of the spectra of the active galaxies^. However, power-law high energy X-rays could 
not be accounted for by this model. Chang and Ostriker^ and Kazanas and Elision^ 
resorted to introducing accretion shocks where the temperature and density would 
be enhanced and the radiation efficiency is also increased. The former introduced 
pre-heating dominated shocks at a very large distance from the black hole, while 
the latter introduced pair-plasma pressure supported shocks close to the black hole. 
Chakrabarti^" ^° showed that centrifugal barrier in a low angular momentum ac- 
cretion flow can produce stable shocks for a wide range of parameter space pro- 
vided the flow has specffic angular momentum everywhere small compared to the 
Keplerian value. This work was further put to test by Chakrabarti and Wiita"'^"'^ 
and Chakrabarti and Titarchuk^^ who showed that shocks could play a major role 
in determining the spectrum of the emitted radiation. Particularly important is 
that, post-shock region which is known as the CENBOL (Centrifugal pressure Sup- 
ported Boundary Layer) and which is the repository of hot electrons, can easily 
inverse Comptonize photons from a Keplerian disk located in the pre-shock region 
and the power-law component of the flow may be formed easily without taking 
resort to any hypothetical electron cloud originally invoked in the literature^'^^^^. 
Thus, the CENBOL region in between the horizon and the shock as introduced 
by Chakrabarti and his co-workers behave like a boundary layer where the flow 
dissipates its gravitational energy. Furthermore, shocks have been found to be os- 
cillating when the cooling is introduced and this is explained to be the cause of the 
quasi-periodic oscillations in X-rays^^. 

However, astrophysical shocks also play a major role in acceleration of par- 
ticles. It is well known that the high energy cosmic rays arc produced by shock 
acceleration^^"^^. These shocks are transient in nature and still play an important 
role in shaping the spectrum. Hence, it is likely that the standing shocks, through 
which majority of the accreting matter must pass before entering into a black hole, 
or forming a jet, should be important to energize electrons. This shocks are very 
stable and remain virtually in place even after non-axisymmetric perturbations^^. 
It is thus fltting that we be interested to understand how the energetics of the high 
energy electrons is affected by these shocks and how the spectrum is affected by 
the synchrotron radiation produced by these energetic particles. 

In Section 2, we present the basic equations and the relevant parameters for the 
problem. We assume the Paczynski-Wiita^^ potential to describe the space-time 
around a non-rotating black hole. We introduce a major emitting component, the 
accretion shocks. In Section 3, we present the relevant heating and cooling processes 
which are taking place in the accretion flow. In Section 4, we discuss how a broad 
band spectra may be obtained. In Section 5, we present a typical spectrum and its 
components. We also give particular emphasis to high energy gamma-ray emissions 
in hard and soft states. Finally, in Section 6, we present concluding remarks. 
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Fig. 1: Cartoon diagram of a quasi-spherical, low-angular momentum accretion around 
a compact object. Dark shaded region indicates a standing shock formed by centrifugal 
force. 



2. Basic Hydrodynamic Equations 

As far as the flow topology is concerned, we assume it to be thin, axi-symmetric, 
and of conical cross-section. The flow is assumed to be from winds of the companion 
star as in the case of a stellar black hole, or supermassive black holes. We do not 
consider the presence of a Keplerian disk, or, if present it is assumed to be very far 
from the region of consideration. The motivation of our present work is to sec what 
the advective disk emits through bremsstrahlung, Comptonization and synchrotron 
emission. In a future work we shall incorporate the Keplerian disk as well. 

Figure 1 gives a cartoon diagram of the flow geometry we consider in this pa- 
per. The low angular momentum flow roughly moves as a freely falling gas in the 
Paczynski-Wiita^^ potential with the velocity profile: 

v(x) = {x - 1)-'/^ (1) 

where x is the radial distance in units of the Schwarzschild radius, Vg = 2GM/c^, c 
is velocity of light, G is the gravitational constant and M is the mass of the central 
black hole. The electron number density, determined from the mass conservation 
law assuming pure hydrogen, is given by, 

where the geometric factor O arises because we are assuming conical flow of solid 
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angle fl (instead of Air valid for a Bondi flow), 

O = 47rcos(e), (3) 

where 9 is the angle made by the surface of the flow with the vertical axis. M is the 
mass accretion rate and rUp is proton mass. In presence of pure hydrogen, Thomson 
scattering will be the most dominating scattering process and the corresponding 
optical depth is given by, 



477, m, Va; — 1 



477 m f 
^^^^ = O ^ 2 V. 



dx, (4) 



where rh is the mass accretion rate in units of Eddington rate Msdci 

M 

fn = — . (5) 

M Edd 

We calculate the magnetic field at a given radial distance from the equipartition 
between the gravitational energy density and the magnetic energy density i. c, 

GMp 

8^ = J^y 

Using (1) and (2), the energy balance equations for protons and electrons can be 
written as, 

where, 7 is 5/3 for non-relativistic electron temperatures (Tg < nieC^/k) and 4/3 
for relativistic electron temperatures (Tg > rricC^/k). k is the Boltzmann constant. 
Since protons are much heavier than the electrons, Tp always remains in the non- 
relativistic domain. T and A contain contributions from all the heating and cooling 
processes respectively. Detailed nature; of these terms will be discussed below. 

It has been shown*^"^" that the shocks in a black hole geometry can typically 
occur at around Xs ~ 5 — 300rg depending on the specific angular momentum A. We 
do not explicitly solve for the shock locations in this paper. Thus, instead of using A 
as a free parameter as in Chakrabarti^" ^"^j wc use Xg to be the free parameter. In a 
transonic solution, the shock strength is also computed from the infiow parameters. 
We use i?, the compression ratio, to be a free parameter as well. 

3. Heating and Cooling Processes 

First of all, we neglect heating due to dissipation as far as the protons are concerned. 
Protons lose energy through Coulomb interaction Agp and inverse bremsstrahlung 
Ajfc. So, 

Kp = Aep + Aifc. (8) 
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Here, the subscript p represents protons. The effect of Aj^ is generally much smaller 
than Aep. Electron-proton coupling supplies energy to the electrons (from protons) 
and is given by, 

Ae, = 1.6 X io-i3 feA/^e^»Ao ^2(y^ _ T,)T-3/^ (9) 

TTip 

where, InAo is the Coulomb logarithm, nip and me are the rest masses of proton 
and electron respectively. Electrons are heated through this Coulomb coupling. 

Te = Aep. (10) 

Subscript e represents electrons. 

Cooling terms for the electrons include bremsstrahlung Ab, cyclo-synchrotron 
Aes and Comptonization A^c of the soft photons due to cyclo-synchrotron radiation. 
For the time being we ignore any Keplerian flow on the equatorial plane which could 
also supply soft photons. The effect of this would be to introduce a bump in the 
soft X-ray and a power-law component due to Comptonization^^. The net cooling 
of the electrons is: 

Ae ^ Ab + Acs + Arnc- (H) 

Explicit expressions for the cooling terms for electrons satisfying Maxwell- 
Boltzmann distribution are: 

Aib = 1.4 X IQ-^'^n^i—Tp) ^'"^ (12a) 
\mp J 

Ab = 1.4 X 10-2^n^Ty2(i + 4.4 x lO-i^Te) (126) 

Ae. = |^fcTe(a;)^, (12c) 
3c^ X 

where Ua is the critical frequency at which the self-absorbed synchrotron radiation 
spectrum is peaked and it can be determined from the relation, 

3 
—I 
2 

where, 

uo = 2.8 X lO^B, (14a) 

de = (146) 

Procedure to determination of Xm is discussed below. 

When the injected electrons obey a power- law distribution, the expressions given 
above will change. For instance, the cooling term due to cyclo-synchrotron photons 
would be given by, 

Ac. = ^aB(-+i)/2(.(?-P)/2 - .(;V/'/^) (15) 



t^a = Tl^oOeXm (13) 



February 2, 2008 15:47 WSPC/INSTRUCTION FILE JD04154 



6 Mandal and Chakrabarti 
where, 

^" mec2(l+p)(3-p)^ 3e ' ' ^ ' 

r(j;/4 + 19/12)r(j;/4 - l/12)r(j;/4 + 5/4) 

^ = r(p/4 + 7/4) • (^^'^ 

and K is the normahzation constant of power-law electron distribution, 

n{£) = K£-P, (17) 

which is obtained using the constraint that the electron number is conserved during 
shock acceleration. The Comptonization is computed by using this cooling term 
augmented by the enhancement factor !F: 

Kmc = Acs^. (18) 

For Comptonization of the thermal seed photons by thermal electrons, T is given 



where. 



m = (20a) 



P=\-exp(-T,,), (206) 

is the probability that an escaping photon is scattered, while, 

A = l+Ae^ + ml, (21) 

is the mean amplification factor in the energy of a scattered photon when the 
scattering electrons have a Maxwellian velocity distribution of temperature 6e, 
rj'2 = 1 — and Xa = hua/mlc^. For Comptonization of the non-thermal seed 
photons (such as generated by the power-law electrons) by thermal electrons, the 
amplification factor can be written as^^. 



where. 



s = ^, (236) 
cl^s = cl^+^, (23c) 
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and 

(?i = ?72 - 1. (23d) 

The amplification factor^^ for the seed photons comptonized by the electrons 
obeying power-law distribution can be written as, 

T = -arRcK^-^ f^, (24) 

3 (3 - p) 

where Rc is the size of the Comptonized region. Smax and Smin are maximum and 
minimum energy of the power-law electrons respectively and ctt is the Thomson 
scattering cross-section. 

In presence of both thermal and noii-th.6rincLl 6l6ctroiis Xui has to be calculated 
by equating combined thermal and non-thermal emission from the CENBOL with 
the appropriate source functions^^ We ignore the effects of power-law electrons on 
the emission of brcmsstrahlung radiation since this radiation is very weak. The 
spectral index is computed by standard procedure^^'^^ . 

4. Solution Procedure 

For a given set of initial parameters, we fix the outer boundary at a large distance 
(say, lO^rg) and supply matter (both electrons and protons) with the same temper- 
ature (say, Tp = Tg = 10^ K). Our result is insensitive to this outer radius as long 
as it is beyond lO^r^ or so. Radial dependence of velocity and density is chosen 
to be those of the freely falling matter (1). Wc then use Runge-Kutta method to 
integrate (7a-b) simultaneously to obtain the electron and proton temperatures as 
a function of radial distance. After we obtained the density and the temperature at 
any point, we compute the radiation emitted by the flow through brcmsstrahlung 
and synchrotron radiation. The degree of interception of these low energy photons 
are computed from the corresponding optical depth these intercepted low energy 
photons are then inverse Comptonized by the hot electrons in the flow. The rest 
are allowed to escape from the flow directly to the observer. We followed the proce- 
dures presented in Chakrabarti & Titarchuk^^ while computing the Comptonized 
spectrum except that our soft-photon source is distributed throughout the flow in 
the form of brcmsstrahlung and synchrotron emission. At the end, we add all the 
contributions to get the net photon emissions from the flow. The geometry of the 
flow is chosen to be conical. The angle Q subtended by the flow surface with the 
z-axis is chosen to be a parameter. 

The shock of compression ratio R causes the formation of power-law electrons 
of slopei9-20: 

p={R + 2)/{R-l). (25a) 

This power-law electrons produce a power-law synchrotron emission with index q 
given by^^ 

g=(l-p)/2. (25b) 
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The power-law electrons have energy minimum at 

^min ~ f^e^^ ^min (26ci) 

and have energy maximum at 



„2 



rripC 



Tmax (266) 



obtained sclf-consistently by conserving the number of power-law electrons and by 
computing the number of scatterings that the electrons undergo inside the disk 
before they escape. This yields, 

r -r + (27) 

max — rmn '■ QR 1/2-' ' \ ' J 

The minimum value of the Lorcntz factor Tjni„ is obtained from the temperature of 
the injected electrons. This temperature is obtained self-consistently through our 
integration procedure. Earlier it has been shown^^ that not all the matter actually 
passes through the shock. Only when the flow passes close to the equatorial plane 
actually passes through a shock. Furthermore, this fraction is time dependent for 
oscillating shocks^*. Thus, in the absence of a fully time dependent solution, we 
assume the percentage of electrons C acquiring a power-law energy distribution to 
be a free parameter. 

We thus have a complete recipe for generating a broad band spectra as a function 
of the following parameters: Q, Xg, R, C and m. The relevant ranges are: < 6 < 90 
5 < a;^ < 300, 0<i?<4, 0<C<1, 10"^ ^ m < 3. The upper limit of R is 
obtained^^ for a gas of polytropic index 5/3: 

Rstrong — T- (28) 

7-1 

In the next Section, we discuss the nature of a typical spectrum obtained for 
a black hole of mass IOMq. Our procedure is equally valid for massive and super- 
massive black holes in Active Galactic Nuclei, and this will be discussed elsewhere. 



5. Results and Interpretations 

5.1. A Typical Broadband Spectrum 

Figure 2 shows the variation of the electron (dotted) and proton (solid) tempera- 
tures (Te and Tp respectively) when rh = 0.1, R = 3.9, 6 = 77°, ( = 0.4 and Xg = 80 
as a function of the radial distance x (measured in units of the Schwarzschild radius 
Vg). Both the axis are in logarithmic scale. It is clear that since Coulomb coupling 
is not very strong, the heating of the electrons is not very efficient. Thus, they start 
becoming cooler closer to the black hole when the number density becomes higher. 
Higher number density increases the cooling of the electrons for x ^ 300. Very close 
to the black hole, especially after the shock at x = Xg = 80, the splitting is dramatic 
and electrons become cooler more rapidly. 
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Fig. 2: Variation of tiie electron (dotted) and proton (solid) temperatures (Te and Tp 
respectively) when rh = 0.1, R = 3.9, © = 77", ^ = 0.4 and = 80 as a function of the 
radial distance x. As flow comes closer to the black hole, the temperatures separate due 
to fast cooling of the electrons due to synchrotron radiation. 

In Figure 3, we present a typical spectrum with all the contributions from the 
accretion flow. The parameters chosen are the same as above. Here, different curves 
are marked with a number. 

The curve marked '1' is the synchrotron emission from the pre-shock part of 
the accretion flow. The curve marked '2' is the Comptonized spectrum of pre- 
shock synchrotron photons due to thermal electrons. The curve marked '3' gives the 
bremsstrahlung emission from the pre-shock region and curve marked '4' is the cor- 
responding Comptonized spectrum. The curve marked '5' indicates the synchrotron 
spectrum from the post-shock region. Comptonized spectra of the synchrotron radi- 
ation from the post-shock accretion flow due to thermal and non-thermal electrons 
are indicated in curves marked '6' and '7' respectively. The curve marked '8' and '9' 
represents respectively the bremsstrahlung emission and it's Comptonization from 
the post-shock flow. The curve marked '10' is the total broad-band spectra from the 
pre-shock and the post-shock regions, fa indicates the synchrotron self-absorption 
frequency. 

It is clear that the power-law electrons, which are the hall-mark of the shock 
in accretion can leave its signature on the emitted spectrum. The curve '6' is soft 
and thermal component, while the curve '7' is hard and is generally from the non- 
thermal power-law component emitted by the power-law electrons. 

When the shock-location is varied, the following changes are expected: when the 
shock is located very close to the black hole, the hot pre-shock flow is the dominant 
source of Comptonized photons while when the shock is very far away, the post- 
shock region (CENBOL) also contributes, but as a whole, the contribution is much 
lower as the outer flow is cooler. When the compression ratio is varied, the optical 
depth in the post-shock region is changed. For a weak shock, the jump in density 
p+ = Rp-, where -|- and — denote the post and pre-shock values respectively, at 
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Fig. 3: A typical broad-band spectrum witii all the contributions from the accretion flow 
is shown for the same parameter as in Fig. 2 (see text for details). 

the shock is not so high. In this case, the power law electrons are not energetic 
enough to leave its signature in the spectrum. When the accretion rate increased, 
it is expected to increase the density of the flow and it becomes difficult to cool 
the matter by the soft photons of the synchrotron radiation. The Comptonized 
spectrum becomes harder. At lower accretion rate, it is easy to cool the flow and 
the spectrum becomes softer. When the percentage of power-law electrons, namely, 
( is increased the spectrum becomes harder at high energy. 

5.2. Spectra in Soft and Hard States 

We now discuss the basic changes in the flow properties which must occur if the 
spectral state is to switch from hard to soft state (in soft-X-ray region) and vice 
versa. In Fig. 4, we present the photon spectra in the hard and the soft states. The 
hard state was created with rh = 1.0, R = 4, Xg = 180, 6 = 80 and ( = 0.01. 
In the soft state, the accretion rate m = 0.1, Xg = 50 and C, = 0.7 was chosen. 
Other parameters were kept fixed. Here too, we did not consider the presence of a 
Keplerian disk, which can also supply soft photons, as we are interested only in the 
behaviour of an advective flow. 

In the soft state, the spectral slope is high and the energy spectral index at 
{F{E) oc i?^"*) from the photons from thermal electrons is found to be at = 1.47 
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Fig. 4: State transitions occurring in a black hole candidate when the accretion rate and 
the shock location is varied. The thermal bump in the hard state causes two intersections 
to form one at a few keV and the other at a few hundred keV (see text for details). 

which that for the not-thermal electrons is a„t = 0.53. In the hard state at = 0.45 
and ant = 1.35. Because of the thermal bump at around lOOkeV, the two spectra 
intersect at two energies, one at lower energy (~ 4— lOkeV) and the other at a high 
energy ~ 300 — 400keV. Normally, the intersection at the lower energy is a pivotal 
point^^. In Fig. 5 we zoom the high energy region of the spectrum and show the 
nature of the each components. The thermal component in the hard state is from 
the thermal electrons in the pre-shock and the post-shock flow. The non-thermal 
power-law component in the hard state is coming from the power-law electrons in 
the post-shock region. The non-thermal component in the soft-state is from non- 
thermal electrons. 

If one compares with observation one notices that the black hole candidates Cyg 
X.130-31 GROJ1719-24 and GROJ4022-h32 all show the similar characteristics as 
presented in Fig. 5 (see, Ling & Wheaton'^^ and references therein). We therefore 
believe successful reproduction of the salient features in our solution shows that 
the non-thermal emission, at least in these class of black holes, is probably coming 
from the accretion flow itself. 

In Fig. 6, we have presented the variation of the spectral indices of the thermal 
component and the non-thermal components as function of the sub-Keplerian rate 
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Fig. 5: The spectra in soft and hard states are shown only at high energy region. The 

thermal bump is due to thermal electrons from the post-shock region while the non- 
thermal power-law components are due to non-thermal electrons generated by the shock- 
acceleration. 

m. The other parameters are, R — 3.9, Xg — 10.0, O = 77, and ( — 0.7. Generally 
speaking, both the spectral indices decrease initially with the accretion rate but 
after about an Eddington rate, where the optical depth itself started getting larger 
than unity, the spectra start becoming softer as expected. 

6. Concluding Remarks 

In this paper, we have explored the way a shock in an accretion flow may be 
identified by observing the spectrum. We have used the shock location and its 
strength as free parameters (at the expense of the specific angular momentum 
and viscosity parameters) although we have used reasonable values discussed in 
earlier analytical studies. We considered the soft photons due to bremsstrahlung 
and synchrotron radiation as the seed photons for the Comptonization. We also 
included the shock acceleration of the electrons and their effect on the emitted 
spectrum. We have incorporated the splitting of temperature of the electrons and 
protons due to radiative processes. 

Our conclusion is that there are several ways a shock may be distinguished in 
the spectrum. At a strong shock, the power-law electrons are produced with a very 
high jmax and that produces a power-law feature in the spectrum. These power-law 
features are seen in both the hard and the soft components in the same manner 
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Fig. 6: Variation of the spectral slopes of the thermal and the non-thermal components 
at and ant respectively as a function of the dimensionless accretion rate. 



as observed in the galactic black hole candidates such as Cyg X-1, GROJ1719- 

24 and GROJ4G22+32. Wc find that since the shocks are inbuilt in the advective 
disk solutions, and electrons can easily produce high energy radiation first through 
shock acceleration and then through synchrotron emission, one need not invoke 
extra source of high energy photons as are usually done in the literature^'^"^*. The 
outflow that is derived from a CENBOL itself can have shocks in jets which in turn 
may also emit power-law emissions and contribute to the total spectra. We shall 
explore this aspect of the problem in near future. 

In Chakrabarti & Titarchuk"'^^ soft photons due to a Kcplcrian disk was respon- 
sible to cool down the post-shock region. In the present paper, we deviate from that 
paper in the sense that we assume that the Keplerian disk is either located very far 
away, or non-existent. The soft photons are locally generated due to thermal and 
magnetic bremsstrahlung processes. In presence of a Kcplcrian disk, the soft X-ray 
bump will be produced at around a few KeV, while for shocks the bump is located 
in UV. Thus, it is possible that these two distinct bumps can be distinguished 
through observation. This work will also be presented elsewhere. 

There has been mush discussions about the so-called advection dominated ac- 
cretion flows or ADAF (see, Esin et al.^^ and references therein) which also claims 
to produce quasi-spherical accretion accretion as the transonic flows described in 
Chakrabarti^^. However, ADAFs or any of its modifications do not have any shocks 
as in our transonic solution and therefore they need to invoke external sources of 
non-thermal electrons. Thus our solutions require lesser parameters that the class 
models based on ADAF. 

This work is partly supported by a RESPOND project funded by Indian Space 
Research Organization (ISRO). 
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